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A perturbative method is developed to calculate the finite size corrections of the low lying 
energies of the asymmetric XXZ hamiltonian near the stochastic line. The crossover from isotropic to 
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The mass gaps of the growth model are found as a function of the growth rate and the substrate 
slope. Higher order corrections to the growth model mass gaps are also calculated to obtain the first 
terms of the KPZ to Edwards- Wilkinson crossover scaling function in the large argument expansion 
in the zero slope sector. 
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I. INTRODUCTION 

Understanding and classifying dynamic scalings of self-affine growing surfaces has been an area of much interest [1]. 
One of the universality classes which attracted much attention is that of Kardar-Parisi-Zhang (KPZ) [2]. The KPZ 
equation is a nonlinear stochastic equation of motion for the height of growing surfaces whose local growth velocity 
depends nonlinearly on the local slope, and is equivalent to the noisy Burgers equation. When the nonlinear term 
is absent, the equation reduces to the diffusion equation termed as the Edwards- Wilkinson (EW) type. There are a 
number of microscopic lattice models which are believed to be described by the KPZ equation in continuum limit, 
e.g. the Eden model, the ballistic deposition model and the restricted solid on solid model [1]. However, identification 
of their universality class has been largely based on qualitative arguments and numerical evidences. An exception 
is the single-step, or the body-centered solid-on-solid model introduced by Plischke et at [3] in one-space one-time 
dimension to which this work is addressed. In this model, the height of each site must differ by ±1 from its two 
neighbors and obey the following dynamic growth rule: Select one of the sites at random; if it is a local valley, i.e. if 
both of its neighboring sites are higher, its height is increased by 2 with the deposition probability pd', if it is a local 
hill, its height is decreased by 2 with the evaporation probability p e ; otherwise, the height does not change. Due to 
the solid-on-solid condition, height configurations of this model have a natural representation in terms of the step 
variable {erf} (af = ±1). If one interpret the sites with a, = 1 as those occupied by a lattice gas particle, the model is 
that of a driven diffusive lattice gas [4]. The time evolution of the probability distribution for the state of the model 
is generated by the A = 1 line, or the stochastic line, of the following asymmetric XXZ hamiltonian [5,6]; 

N ~ 

n = - «^ - ^yV^r+i - ^^r^+i} (i) 

8 = 1 

where N is the number of sites, af are the spin 1/2 raising and lowering operators, Cjy- +1 = a" and s = (pd~Pe)/(Pd + 
p e ) is the growth rate, m = X^fci a i I ^ is the mean slope of the substrate which is conserved in this model. 

The hamiltonian with A as a parameter has other physical applications as well. For example, it describes endpoints 
of facet-ridges in equilibrium crystal shapes and persistent current transitions in mesoscopic metallic rings [6-8]. The 
hamiltonian is an anisotropic limit of the transfer matrix of the general asymmetric six vertex model. When s = 1, 
the latter becomes the five- vertex model which describes, among others, an interacting domain wall system [9]. 

When A < 1, the hamiltonian is critical and possesses the conformal invariance [10]. Thus the mass gaps, which are 
the differences of low lying energy eigenvalues from the ground state energy, scale as l/N , the proportionality constant 
being the conformal dimension multiplied by a non-universal anisotropy factor. At the stochastic line, the ground 
state of the hamiltonian becomes trivial with zero energy in each sector, a subspace of % with fixed value of m. The 
mass gaps in this case are expected to scale as l/N 3 ! 2 , since the dynamic critical exponent z of the KPZ universality 
class are known to be 3/2 in one dimension [11]. That this is indeed so has been demonstrated by Gwa and Spohn 
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[5,12] from an analysis of the Bethe ansatz equation for the special case of s = 1, and m = 0. Also, recently Neergaard 
and den Nijs [6] have analyzed numerically obtained mass gaps to investigate the crossover behaviors from the KPZ 
line < s < 1 to the s = EW point where the mass gaps scale as l/N 2 . The bulk properties of the asymmetric 
six-vertex model at and near A = 1 are studied by Bukman and Shore [13,7] in connection with its application to the 
equilibrium crystal shape. Correlations in the asymmetric XXZ chain with reflecting boundary conditions which has 
finite energy gaps in the bulk is studied by Sandow and Schiitz [14] and by Henkel and Schiitz [15] in the context of 
driven diffusive lattice gas system. Also, the continuum limit of Eq. (1) at A = 1 is derived by Fogedby et al. [16] 
from which the dispersion relation u> oc k 3 l 2 is obtained for s > in the quasi-classical limit. However, it is of interest 
to extend our exact knowledge of this model which can play a central role in the non-equilibrium growth problems. 

In this work, we analyze the finite size corrections of the energy E^(A, s, m) of Eq. (1) for low lying levels using 
the Bethe ansatz equation and develop a systematic expansion of En in powers of l/N 1 ' 2 for general s and m and 
for small 1 — A. This method enables one to obtain general expressions for mass gaps associated with various energy 
levels generalizing the work of Gwa and Spohn considerably. In the next section, we review the Bethe ansatz equation 
and put it in the form suitable for analysis for small 1 — A. In Section III, finite sums are expressed as an infinite 
series in l/N 1 ! 2 assuming certain analyticity property of the phase function introduced in Section II. The sum formula 
derived in Section III is then used in Section IV to determine the phase function self-consistently and to express the 
energy of Eq. (1) in the form of a perturbation series in l/N 1 ! 2 . In the next section, we solve explicitly the general 
expressions of Section IV to obtain leading order behaviors of En ■ To the leading order, mass gaps are then obtained 
in terms of universal scaling functions which describe the crossover from the isotropic (A < 1, z=l) to anisotropic 
(A=l, 2=3/2) scaling behaviors. The scaling variable is found to be s _1 (l — m 2 ) 1/,2 (l — A)N 3 / 2 up to an arbitrary 
numerical factor. In particular, the real part of the mass gaps are proportional to s(l — m 2 ) 1 ! 2 N~ 3 I 2 multiplied by 
the universal scaling function. Higher order terms are also determined to 0(l/N 4 ) in the m = sector at A = 1. The 
latter series determines the first three terms of the KPZ to EW crossover scaling function for large argument where 
the scaling variable is sN 1 ! 2 . We then conclude with summary and discussion in Section VI. 



II. BETHE ANSATZ EQUATION 

The asymmetric XXZ hamiltonian Eq. (1) is diagonalized by the Bethe ansatz. This is most conveniently seen from 
the fact that it is an anisotropic limit of the transfer matrix of the general asymmetric six-vertex model. This we 
show in Appendix A for self-containedness. Since Eq. (1) conserves m, we can consider each sector of fixed number 
of down-steps separately. We use the notation Q to denote the number of down-steps (down-arrows in the language 
of the six vertex model) and let 

q = Q/N=(l-m)/2. (2) 
The energy of Eq. (1) in the sector Q is given by 

Q 



■ s -i 1 — s 

j=i 

where the complex fugacities {zj} satisfy the Bethe ansatz equation 

Q 



E N = "£{A-^z^-^-l z . } (3) 



( iW -1 TT 1 + 6 ZiZj 2Ae Zi (a\ 

{ ' J-J- 1 + e- 4H z iZi - 2Ae- 2H Zi ' { ' 



The notation z used in this section and Appendix A refers to the fugacity and should not be confused with the 
dynamic exponent. Different solutions of Eq. (4) give the energy of different eigenstates. Here, we have used the 
standard six-vertex model notations of [17] for H and A which are related to s and A of Eqs. (1) and (3) by 

A = A/cosh(2 J ff) (5a) 

s = tanh(2if). (5b) 

Since Eq. (1) with s < is related to that with s > by the simple particle-hole symmetry, we only consider the 
region < s < 1. We are interested in the region near A = 1 and parametrize A by 
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A = cosh v (6) 
(y > 0). The fugacity variable z is transformed to a more convenient variable x defined through 

z = exp(2ff -v) l ~ X (7) 
1 — ax 

where 

a = exp(-2i/). (8) 

We note here that the five- vertex limit is achieved by letting v and H approach infinity (a — ► 0), keeping 2H — v 
fixed. 

Using these notations, Eqs. (3) and (4) then take the simpler forms, 



E 



sinh v v^-\ , Xi ctXi 



N 



cosh 2H 

.7=1 



1 _ \ W Q _ 

2H-V 1 ^ I _ ^ 1)0-1 ^ aX j ^Q-j 



respectively. Next we introduce, following de Vega and Woynarovitch [18], the function of complex variable 

iZ N (x) = (2H + log }~ X + 1 f^logz, - log l - aX3 ' X ) (11) 

xi(l — ax) ^ ~[ 1 — ax/Xj 

and its derivative 

R N (x) = ixZ' N (x). (12) 

Zn(x) is called the phase function in this work and plays a central role below. It depends on {xj} , H , v and q. Eq. 
(10) then can be written as 

Z N(x j )= 2 ^rIj (13) 

where Ij are half-integers (integers) for Q even (odd). The set {Ij} specifies a particular energy level and is related 
to the momentum of the state by 

j = l 3=1 

It is well established that the ground state is obtained if {Ij} is chosen as 

I i = - < ^p L +3 for i = !.•••. q ( 15 ) 

while low lying excited levels are obtained by choosing sets of {Ij} differing from that of the ground state in a few 
places near the end. For example, at A = 1, the first excited state is associated with the set 

Ij = ~ ^ ^ 1 + J for i = l,--,Q— 1, and Iq = ® ^ 1 (16) 
and has momentum P = 2tt / N , while the first excited state in the P = subspace is associated with the set 

7i = — , Ij = \-j for j =2,...,Q- 1, and J Q = — — . (17) 
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In the thermodynamic limit N —> oo, the roots xj for the ground state and its neighboring states form a continuous 
contour in the complex- a; plane. It is the inverse image, Z^ 1 ((/>), of the segment of the real axis — irq < <f> < irq of 
the complex <j> plane. Here Z^ 1 denotes the inverse function of Z^ . R 00 (x)/2ir is then the density of roots along 
the contour. In particular when A = 1, the contour turns out to be a closed contour starting from e™g/(l — q) and 
ending at 

x° c =e-"q/(l-q) (18) 
enclosing the origin clockwise, as <f> varies from — irq to irq [13,19]. In other words, we have 

Z oo (x c ) = TTq and Z x (>°e 2 ™) = -Trg. (19) 
The explicit form of Z^ for A = 1 is 

iZ^x) = log 1^ + log(g*(l - g) 1 "'). (20) 

This is obtained by integrating the result for R^^x) of [13] which is R^^x) = —q — x/(l — x) in our notation and 
by fixing the integration constant using Eq. (11). The branch cut for the log function is along the negative real axis. 
An actual form of the contour is shown in Fig. 1 for the case of q = 1/2. It is to be noted here that it has a cusp at 
x = x° c due to the fact that Z'^^x®) = 0. Therefore, if one performs a power series expansion of Eq. (20) at x° c) the 
first order term vanishes and the result becomes 

Z^x) = ±Ttq - i (1 ~ q) \ x - x° c f + 0((x - x° c f) (21) 
Iq 

where the +(— ) sign is for x below (above) the branch cut. Accordingly, Z^ 1 ^) has the square-root type singularity 
at <f> = ±7rg and takes the form 



Z~ 1 (Trg -Z) = x °-J ^h-^ y/i£ + 0(Q (22) 



near the upper end of <f> and 



z^(-^q + = ^ + ^77^^ +0(0 (23) 



(l-g) 3 

near the lower end of <f>. We shall see later that this square-root singularity is the origin of the singularity of free 
energy as a function of A [13] and the unusual 1/ 'TV 3 / 2 scaling of the mass gaps [5]. 



III. SUMMATION FORMULA 



In this section we assume that the phase function Z^(x) is known and derive a general expression of finite sums 
in terms of its expansion coefficients at a critical point x c . The problem at hand is then to evaluate sums of the form 
(See Eqs. (9) and (11).) 

Q 

S[f]=J2f( x s) 

i=i 

= E/(^ 1 (^-)) 

i=i 

Q 2?r 1 

= J2f( Z N 1 (-^ + 1 Tf(J- g))) + other terms (24) 
i=i 

where ' other terms' above denotes the terms which differ from that of the ground state; for example, for the first 
excited state characterized by Eq. (16), 
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other terms = f(Z^\*q + £)) - f(Z^\*q - ^)), (25) 

while more terms are added and subtracted for other levels. In the following, we assume S[f] is for the first excited 
state for the sake of notational simplicity and comment on other levels when appropriate. 

To evaluate the finite sum in Eq. (24), we use the sum formula employed by Gwa and Spohn [5,20]; 

E F(J) = [ Q F(t) dt + \{F{Q) + F(l)) + 2 F{Q ^l;^ ht) dt (26) 



l ^ Jo 



where F(s,t) = (F(s + it) — F(s — it))/2i. The finite size corrections of the sum are determined from properties of 
the summand near the end points of the sum. To produce the l/N 3 ! 2 scaling correctly, the effect of the square-root 
singularity discussed above should be taken account. Thus, we first assume that Z^(x) has a vanishing derivative at 
a certain value of x and define x c by the relation 

Z' N (x c ) = 0. (27) 

For large N, x c is expected to be close to x° c but may pick up an imaginary part in which case we let the branch 
cut of log function pass through x c . Next we assume that Z' N (x) is analytic at x c . This assumption is based on the 
observation that R^(x) = — q — x/(l — x) exactly in the five-vertex limit and that expansion of Eq. (11) in the power 
series of a does not introduce extra singularities. Also define 6± by the relation 

Zn(xc) = irq — ib + and Z^(x c e 27 ' 1 ) = — 7rg — i6-. (28) 

Comparing Eq. (28) with Eq. (19), we also expect 6± to be small for small 1 — A and large N . Presence of such 
terms is crucial in the expansion method we are describing below and has been anticipated from the solutions of the 
five- vertex model by a different method [21]. 

Definitions of x c and 6± together with the assumption of analyticity of Z' N (x) at x c then allow us to put Z^ 1 in 
the form 



ZK 1 (irq-0 = Xc+Y j a m (-iVh+^r (29) 



for the upper end and 



Z- N \-^q + = x c + J2 a m (iy/6- - (30) 

m = l 

for the lower end of the sums in Eq. (24). Eqs. (29) and (30) implies that Z^(x) has the expansion 

iZ N (x) = ±irqi + 6± + -^(x - x c ) 2 - ~^-( x ~ x cf + ••• ■ (31) 

x c , 6± and a m are to be determined self-consistently. They are in general complex but are defined in such a way that 
their leading term in the following perturbation expansion is real. Before proceeding further, we argue a posteriori that 
=8-. The following steps of the derivation of Z^(x) show that it is of the form iqlogx + (parts analytic at x c ) 
regardless of 8+ = 6- or not. Then Zn(x c ) — Z-^(x c e 11[% ) = 27rg implying 5+ = 6- from Eq. (28). Hence, we put 



(32) 



from now on. 

Then, the summand of Eq. (24) near the upper and lower ends can be written in the form 



f(Z~ 1 (vrg - 0) = f(x e ) +J2 A ™ W H V6 + %r - ( 33 ) 

m = l 



fiZ-'i-Trq +0) = f(x c e 2 ™) + ]T A m [f] (i^3^)™, (34) 
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respectively. In the above, the coefficient A m [/] is given by 



m 7, 

A m [f] = J2-^f (k H*c) (35) 



k = l 



where ^ (1 < k < m) is related to a m through the relation 

oo oo 

(J2 a ^ m ) k = J2 b ™- kem > ( 36 ) 

m = l m — k 

and /( & ) denotes the fc-th derivative of f(x). When Eq. (26) is applied to Eq. (24), the first term on the right hand 
side of Eq. (26) contributes 

— / KZ- 1 (</,)) dt = -— f f(x)Z' N (x) dx 

+ ^{ /(^W)^+ /(^W)^ (37) 

L ~ yJirq-ii J-wq + Tr/N J 

where § denotes the integral along the closed contour in a;-plane in the counterclockwise-sense. Using Eqs. (33) and 
(34) for the correction terms in Eq. (37) and the rest of the terms in Eq. (26), one can easily see that the m-th order 
terms in the series of Eqs. (33) and (34) contribute to order N~ m l 2 in S[f], provided the variable 

7T 

is fixed to a finite value. To collect the 0{N~ m l 2 ) contributions, we are thus led to define a set of scaling functions 
Y m (y) as follows; 



Y m (y) = Y^(y) + (-i^i)™ - (-i^Ti)-, (39) 

where 



Y m (y) = Re^— TY (- i Vy + i ) +,J o dt) (40) 

is the contribution from Eq. (26) and the last two terms in Eq. (39) come from the two terms of Eq. (25). Useful 
properties of Y m and are discussed in Appendix B. In terms of Y m 's thus defined, Eq. (24) then can be written in 
the form of a series 

S[f] = lf(x)R N (x) - + Uf(x c ) - f(x c e 2 ")y+ £ A m [f] Y m (y) e m (41) 

/7n / xl z — ' 

J m = l 



where e is our perturbation expansion variable 



and i?jv is given by Eq. (12). The second term in Eq. (41) gives a contribution Try for f(x) = logx but is zero for / 
analytic at x c . 

IV. GENERAL FORM OF Z N AND E N 

The sum formula Eq. (41) is derived pretending that Zn{x) is known. One then applies Eq. (41) to Eq. (11) which 
is the defining equation of Z^ to determine it self-consistently. Here, f = f z where 

/zOO = logs' -log (43) 
1 — ax/x' 
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This then leads to an integral equation for R^(x) which is solved by Fourier method. To handle the contour integra- 
tions, we represent fz as a Laurent series 



o-H r' 

Mx') = \ogx' + J2 (-)"• (44) 



This is justified since < a < 1 and integrations can be made keeping \x' / x\ = 1. Using Eqs. (41) and (44) in Eq. 
(11) then gives us 

1 f dx' 
iZ N (x) = -q\ogx + G(x)- — i f z (x')R N (x') — (45) 

Ziri I x 



where 



G(x) = (2H -v) + log i + 6 

1 — ax 



m = l fc = l 



Here and below stands for the sum over all integers n. Let us denote the Fourier mode of G(x) by G n ; 

G(x) = Y J G n x n . (47) 

n 

Taking derivative of both sides of Eq. (45) with respect to x, one then obtains the integral equation for Rjy(x); 

R N (x) = -q + V nG n x n - -L I f R ( x <) R N ( X <) ^ (48) 
A — ' Ziri I x 

n 

where 



fR(x') = -J2 aH (-T- (49) 



The solution of Eq. (48) is 



Putting this back into Eq. (45), one obtains 



R N{x) = - q + Y J ^ R x\ (50) 



iZ N (x) = wqi -qlog— + G + J2 tS\\ (*" ~ <)■ ( 51 ) 



Inserting explicit forms of G n , we then finally obtain 



oo m (—\) k 

iZ N (x) = wqi + (2H - v) + 6 - - V V b m>k y ~ lr J k -Y m {y) e m+2 

7T fC X . 

m = l & = 1 

oo m 



g log — + log - 1 — — + - E E W {fffcfa) " ^K)} Vm(!/) £ m+2 (52) 



X Q 1 eZ? Q '7T 

m = l fc = l 



with 



a \n\ 

= J k £(" + + 2) • • • (n + A - 1)- R (-)". (53) 

k\x"i z — ' 1 — a> n > x r 

The series in Eq. (53) is convergent for a < \x/x c \ < l/a. In deriving Eq. (52), we assume \x\ < 1 which is valid 
for q < 1/2. However, final results can be extended to the region q > 1/2. Z^ in the five- vertex limit (a — ► 0, 
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2H — v=— log A) differs from that of the general six- vertex case only by the last double sum of Eq. (52) which is 
To determine 6 (i.e. y), x c , and a m 's, we use their defining equations in Eq. (52). Eq. (28) then gives 



2H - V =* EEW^™r +2 , (54a) 

m=l k=l c 

Next, comparing iZ^\x c ) derived from Eq. (52) with the coefficient of (x — x c ) k of Eq. (31), we obtain an infinite 
hierarchy of equations; 



oo m 



= - ^ - ~, ~ +-EEWJ1 ( x c)Y m (y)e m + 2 , (54b) 



q 1 



m = l fc = l 



oo m 



a;? (1 — I,) 2 7r ^-^ ^-^ 

c v c/ m = l fc = l 



-12a 2 -2g 



1 oo m 

- EE^^n^^tr 12 , (54d) 



af a;? (1 — x r ) 3 it 

1 c \ c; m = 1 k = 1 

etc. A merit of Eq. (54) is that it can be used to determine y, x c and the a m 's order by order in e. 
Having determined, though formally, Zn{x), we next apply Eq. (41) to Eq. (9) using f = f E where 

oo 

/^ = "E( 1 (55) 

n = l 

If Rn(x) is of the form of Eq. (50), then a simple application of the residue theorem gives 

If dx °° 

— f f E (x)R N (x) — = y~] nG- n . (56) 
Ziri I x z — ' 

n = l 

This simple result is due to the particular form of j E . Combining this with the explicit form of fg\x c ), we find 

m=l k=l v ' 

It is convenient to eliminate the first term in the double series of Eq. (57) using Eq. (54a) to make the series start 
with 0(s 2 ) term. Thus we write En as 

^-^{ra^-*"-^} (58a) 

where 

oo m 

T * = J2z2 Ck [* J w Y ™ (^) £m ( 58b ) 

m = 2 k = 2 

with 

C k [x e ] = ( c ' • (58c) 

L J jfe(l - xJWx*- 1 K ' 

Eq. (58) together with (54) is the central result of this work. For energy levels other than the first excited one, one 
needs to replace Y m (y)'s in Eqs. (58) and (54) by appropriate ones; e.g. by Y^(y) for the ground state energy. Y m 's 
for the energy level of Eq. (17) will be discussed at the end of the next section and in Appendix B. 
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V. LEADING ORDER BEHAVIORS OF THE MASS GAPS 



We now find perturbative solutions to Eq. (54) and discuss their consequences. To the zeroth order in e, x c and a m 
take the values x° c and a^, respectively, where x° c is given by Eq. (18) and a^'s are obtained by solving Eqs. (54c), 
(54d), etc. with s = or, more conveniently, by inverting the Taylor expansion of Z^x) around x° c) Eq. (21). The 
first few a® 's are 



_ / 2g 2(1 + 3) _ l + llg + ^ 

a i-y(T^)I' 2 ~ ~3(T^qf and a3 " ai 18 g (l- g ) • (59) 

Inspection of Eqs. (54c), (54d), etc. shows that a m = + 0(s 3 ). Furthermore, the 0(s 3 ) term of Eq. (54b) 
is proportional to g'i(x c ) which can be shown to be identically zero from its definition, Eq. (53). Thus, we have 
x c = x° c + 0(e 4 ). Therefore, using x c = x° c and a m = a° m in Eq. (54a) and (57), we get correct expressions to 0(s 5 ) 
and 0(e 3 ), respectively. Also, since 2H — v = 0(e 3 ), and A = 1 — s(2H — v) + 0((2H — v) 2 ) from Eq. (5), we can 
put 2H — v = (1 — A)/s correct to order s 5 . Taking these considerations into account and using Eq. (36) , Eq. (54a) 
then becomes 



1 - A = I , / — ^ Yi{y) £ 3 + -L^l-Y 2 s* + J 2 1 q + q ] Y 3 (y) e 5 + 0(e e ) (60) 

^ {]] (li 1 - v 33(1-3) y 3(1 - 3) 183(1 - 3) j 

which can be used to determine y as a function of (1 — A). This suggest that the appropriate scaling variable is 



To the leading order in e , y is related to t by 

t = Y^y). (62) 

Similarly, Eq. (58a) becomes 



E N = -3(1 - 3)(1 - A)N + s{(l - 23) Y 2 e 2 + -^/2q{\ - q)Y 3 {y) e 3 + 0(e 4 )}. (63) 

which gives En as a function of A, s = tanh2if and 3 to leading orders. Here, Y2 = 2i is related to the momentum 
of the level by Eq. (B7) and y is a function of the scaling variable t through Eq. (62). 

The mass gap is the difference of En from the ground state energy E° N . Since the level dependences of En enter 
only through definitions of Y m (y), E° N is obtained by replacing Y m by Y^ in Eqs. (63) and (62). Thus we have 

E N -E° N = s(l - 2q)27ri/N + 2 -{Y 3 (y) - Y 3 °(y') }s^2q(l - 3) (^f 2 + 0(l/N 2 ) (64) 

where y is the solution of Eq. (62), y' is the solution oft = Y°(y'), and Y 2 ° = is used. We remark that Y^y) — Y^(y') 
is a function oft only with no other parameters. Thus we define a universal scaling function as 

Ht) = ^ 3/2 {Y 3 (y)-Y 3 (y')} 

= ^ 3/2 { W^)) " Y^Y? "!(*))} (65) 

to put Eq. (64) in the form 



E N -E% = 2Trism/N + s\/l - m 2 f(t)/N 3/2 + 0(1/N 2 ) (66) 

where we recall from Eq. (2) that m = 1 — 2q. is the universal scaling function describing the crossover from the 

KPZ 1/N 3 ! 2 scaling to the isotropic 1/N scaling as will be shown shortly. It is universal in the sense that the real part 
of the mass gap, scaled by the factor s\/l — m 2 / 'TV 3 / 2 depends on s, m, 1 — A and N only through the combination 
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given by Eq. (61) as N —> oo. T(i) is shown in Fig. 2 for < t < 4. It is not defined for t < 0. Higher order terms in 
Eq. (66) give the correction to scaling. 

For the single-step model, we have A = 1, which implies t = 0. In this case, E° N is exactly since % generates 
a stochastic process. This is reflected in Eq. (65) by the fact that when t — ► + , y' — ► — oo and Y 3 °(y') — ► 0~ as 
discussed in Appendix B. Thus we have JF(0) = y27r 3 / 2 Y 3 (t/ )/3 where y is the solution of Yi(y ) = 0. Numerically, 
T(Q) = 6.50918933794 . . .. This is exactly the amplitude of the mass gap obtained by Gwa and Spohn [5] for s = 1 
and m = 0. The dynamic exponent z is then 3/2 for any s and m. Further implicatons of Eq. (66) are discussed in 
the next section. 

When 1 — A is finite and positive, t — ► oo for large N . In this case, the model is in the critical phase described by 
conformal field theory with central charge c = 1 and hence is expected to possess mass gaps which scale as l/N [22]. 
To check this, we use the large y behaviors of Y m (y) and Y^(y') discussed in Appendix B to obtain Y 3 (y) and Y 3 °(y') 
as a function oft. The result is 

2 /3A 5/3 1 /3A 1/3 



yM) = - 5 {j) -e(T) +0(r ' , • (67a) 

*<« =-i (?)"'+ t (If + -'- 1 '- 



From Eq. (67), we have 



^^(t)"W.> 



(68) 



as t — ► oo. The first terms in Eq. (67) contribute to the O(N) bulk term in Eq. (63), while the second terms give the 
0(1/ N) correction terms. Using Eqs. (68), (61) and (2) in Eq. (66), we then have 

E N -E° N = ^ |«mi + (^s 2 (l - m 2 ) 2 (l - A)) {1 + 0((1 - A) 2 / 3 )} + 0(1/ N 2 ). (69) 



Also, using Eqs. (67), (61) and (2) in Eq. (63), we have the bulk energy 

~ 2/3 

El/N = --(1 - m 2 )(l - A) - 4z f — ' 

JV^oo 



v Iim E%/N = -1(1 - m 2 )(l - A) - 1 ^ ^ S ~^(l - m 2 ) 4 / 3 (l - A) 5 / 3 {1 + 0((1 - A) 2 / 3 )}. (70) 



The 0((1 - A) 2 / 3 ) corrections in Eqs. (69) and (70) come from higher order terms of e not shown in Eq. (63). We 
note here that both the mass gap scaling and the ground state energy singularity in 1 — A is determined by the 
e 3 ~ N~ z term of Eq. (63). The t 5 / 3 behavior of Y^(y') and Y^(y) in Eq. (67) is necessary to produce O(N) terms 
when combined with the s 3 factor. Therefore the ground state energy singularity exponent 5/3 is related to the 
dynamic exponent by 1 + 1/z as noted in [7] for the six vertex model. We also note that these fractional powers 
result from the fact that e ~ TV -1 / 2 which in turn originates from the square-root singularity of Z^//(<f)) at <j> = ±irq 
as mentioned in Section II. 

To the order discussed so far, the growth rate (s) dependence of the mass gap is rather trivial. Recently, Neergaard 
and den Nijs [6] have studied the mass gaps in the q = 1/2 (i.e. m = 0) sector numerically to discuss the crossover 
from the KPZ l/N 3 ! 2 scaling to the EW l/N 2 scaling as s — ► for A = 1. To obtain nontrivial s-dependence, we need 
to include higher order corrections in Eq. (58). As discussed above, x c and a m differ from their zeroth order value 
by 0(s 4 ) and 0(s 3 ), respectively. We use these zeroth order values in the double sums of Eq. (54). This produces 
errors of 0(e 7 ) in Eq. (54b) and 0(s 6 ) in Eq. (54c), (54d), etc. We then solve for x c and a m perturbatively from 
the approximated Eq. (54) to obtain x c correct to 0(s 6 ) and a m to 0(s 5 ). Using those values in Eq. (58b) would 
give correct result to 0(e 7 ). However, there occurs an extra simplification for q = 1/2. The first term of Eq. (58b) 
is proportional to C2[« c ]ai which becomes simply 1 — 2q at the zeroth order. Thus the knowledge of a\ to 0(s 5 ) is 
sufficient to determine it to 0(s 6 ), resulting in Tjy correct to 0(s 8 ). Symbolic manipulation softwares are used for 
the algebra. The result for Tjy at q = 1/2 is 

Tn = ^Y 3 (y)s 3 - ^§Y 5 (y)s 5 + ^ 2 ( a )(Y 1 (y)Y 3 (y) - Y 2 )s 6 + ^Y 7 (y)s 7 

+ -{Ufo(<*) + 8/?4(a))Y 3 (j/) 2 - 16Ma)Y 2 Y 4 (y)}s 8 + 0(s 9 ) (71) 

7T 6 
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where (k even) are functions of v defined by 

_ _ r)k r\/\ n \ 

fc( tt ) = £r^R ( 72 ) 

n 

(See Eq. (8).) and y is to be determined from Eq. (54a). When A = 1, the equation for y to the necessary order for 
q = 1/2 becomes 

= Yi(!/) + \Y 3 (y)e 2 + ~^My)e A + 0(s 5 ). (73) 

The solution of Eq. (73) can be obtained in the form y = y + yie 2 + y 2 s 4 • • • . Using it in Eq. (71) and using the 
relation En = sTjy for A = 1 , we then finally obtain 

= - - +{i><»>+ 7 4£$ « 7 

1 39 

+ -{^y 3 (j/o) 2 - WY 2 Y 4 (y )}[] 4 (a)s 8 + 0(e 9 ). (74) 

In deriving Eq. (74), we have used Eq. (Bl) repeatedly. That 0(e 2 ) and 0(e 4 ) terms are missing is a simplifying 
feature of the q = 1/2 sector. This is valid for A = 1 only and fa's are now functions of s since 

a = exp(-4 J ff) = (75) 
when A = 1. Thus the mass gap for A = 1 and q = 1/2 is of the form 

where and 5j are constants which are found numerically, except Bo- Ao is nothing but ^"(0) introduced above. 
Numerical values of the coefficients are as follows; 

A = 6.50918933794..., A l = 10.8607132611 ... , A 2 = -10.3825789166 ... , 

5 = 32tt 2 , B l = -2407.62590782.... 

This expression for energy is in complete agreement with existing numerical values of En obtained from diagonalization 
of Eq. (1) for small TV's (N < 18) in [6]. For example, A 1 /(A tt) is estimated in [6] be to 0.531±0.002 while it is 
0.531106443911. . . in our result. Moreover, the numerical sequence of (E N N 3 l 2 /s - A - Ai/N)N 2 for #=6,8,. . .,18, 
obtained using the s = 1 data of [6] together with the exact values of Ao and Ai given above extrapolates to around 
— 10.35 ~ —10.40 consistent with the exact value of A 2 . A similar procedure using the s = 0.999 data could produce 
Bo but with less precision [23]. 

It is of interest to consider the KPZ to EW crossover scaling function near s = 0. When s — ► 0, it can be shown 
that /3k defined in Eq. (72) and (75) behaves as 

(3 k (a) -> 2k\((k + 1)(1 - a)- (&+1) ~ kV2~ k C{k + l) S - (k+ ^ (77) 

where £ is the Riemann zeta function. One notes from this that j3 2 {a)/N 3 and ^(^/N 4 in Eq. (76) are both 
0(1 /N 3 ! 2 ) provided s\/ r N is kept constant as s — ► 0. Therefore the natural crossover scaling variable is u = s\/ r N and 
the mass gap takes the scaling form 



E N = g(sVN)/N z (78) 

as s — ► with sy/~N fixed. Eq. (76) then gives the first three terms of the large argument expansion of the crossover 
scaling function Q(u) as 

g(„) = AoM + 4 + 4 + --- ( 79 ) 

u u 
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where B' = C(3)5 /2 and B[ = 3((5)5i/2 . At the EW point (s = 0), it is known that E N = 2ir 2 /N 2 exactly [6] 
and hence we have Q(0) = 2ir 2 . The scaling form Eq. (78) has been anticipated in [6]. However, the actual form of Q 
proposed in [6] is of the form 27r 2 ^/l + u 2 /ir 2 = 2ir\u\ + 7t 3 /|m| + • • • which we find is valid only approximately. 

So far, only the first excited state characterized by the set Eq. (16) is discussed. However, the above results are 
readily applied to other energy levels provided the scaling functions Y m (y) are substituted by appropriate ones. For 
example, to consider the momentum level mentioned in Eq. (17), the Y m need to be changed to 

Y m (y) - Y°(y) + {-iyfy^iT + (i\/y + i) m " " {iy/y^T (80) 

where Y^(y) are still given by Eq. (40) and the remaining 4 terms comes from two excitations and two holes in Eq. 
(17). Here, Y 2 = from Eq. (B7). For q = 1/2, only the numerical values of the Y m (yo) in Eq. (74) change. See 
Appendix B for the value of t/o- The coefficients of Eq. (76) are then found to be 

A = 16.0176269046..., A l = 44.0512971116 ... , A 2 = -7.2545175082 ... , 

5 = 0, B l = 12315.0898394.... 

In [6], Aq/t: 3 ! 2 for this level is estimated to be 2.87±0.01 which is consistent with our exact value 2.87655951907. . .. 
However, their conjectured exact value of 2.8633717. . . is off by 0.5%. Note that the non-trivial s dependence appears 
only at 0(l/N 4 ). Other higher levels can be treated similarly. 

VI. SUMMARY AND DISCUSSION 

In this work, we have developed a perturbative scheme which allows one to calculate the finite size corrections of 
the low lying energies of the asymmetric XXZ hamiltonian Eq. (1) near the stochastic line and obtained their leading 
order behaviors for arbitrary s and arbitrary m. Formally, the energy is given by Eq. (58) together with Eq. (54) 
which determines perturbatively y, x c and a m 's entering Eq. (58). We remind the reader that q is trivially related to 
the substrate slope by Eq. (2) and other notations in Eq. (58) are defined in Section III. 

The mass gap of the first excited level is found to take the scaling form given by Eq. (66) where the scaling variable 
t is given by Eq. (61) and the universal scaling function is defined in Eq. (65) and shown in Fig. 2. behaves 
as ~ t 1 ! 3 for t —> oo and jF(0)=constant. The result can also be applied to other levels. Thus if we use superscript I 
to distinguish various low lying energy levels, we then have the scaling form for the mass gap as 

E l N -E% = ismPt + s\/l ~ m 2 F i (t)N- 3/2 + 0(#" 2 ) (81) 

with Pi=(2tt / N) xinteger the momentum of the level I, Tt the level dependent scaling function. Ti\t) is also defined 
by Eq. (65) but with appropriate Y m 's, e.g. Eq. (80) for the lowest level with PI = 0. This is the first main result of 
this paper. 

When A;$l,2— >ooas./V— >oo and the mass gaps scale with l/N as shown in Eq. (69) for the first excited level. 
In this case, the continuum theory of Eq. (1) is described by the central charge c = 1 conformal field theory. The 
first excited level characterized by Eq. (16) is the first descendant of the identity operator in the sense of conformal 
field theory. Thus the mass gap is expected to behave as 2tt(/N, ( being the complex anisotropy factor which is 
proportional to the modular ratio of the torus on which corresponding field theory is defined [24]. Eq. (69) then 
determines ( for A close to 1, with the imaginary part accounting for the deformation of the square lattice to a 
parallelogram. The ground state energy itself shows a sing ularity of the type (1 - A) 1+1 ^ with z = 3/2. The 
crossover from the conformally invariant region A < 1 to the KPZ line is described by the scaling functions Ti\t) as 
functions of the scaling variable t. 

For the single-step growth model, t = 0, E° N = and Ti{i) becomes a universal constant Ti(0). A quantity of 
interest here is the stationary slope-slope correlation function given by [5] 

G(r,T)=<0\a z o e-^ n a z r \0> -m 2 (82) 

where |0 > is the ground state in sector m, r stands for time and G should not be confused with the notation 
used in Section IV. Our result then implies that the correlation function takes the form, after taking the spectral 
decomposition and using Eq. (81), 

G(r,r) = < 0\a z \£ >< £\a z \0 > exp{-iP^(r + smr) - s\/l ~ m 2 ^(0)|r|/# 3/2 } (83) 
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as r — ► oo, N -> oo with t/N 3 / 2 fixed. H ere, the sum is over excited levels whose mass gap scales as l/N ' , 
\l > stands for excited state and the relation <r z = exp(— iVr)<j z Q exp(i'Pr), V being the momentum operator whose 
eigenvalue is Pi, has been used. Therefore decay of the correlation is governed by a characteristic time which scales 
as N z with z = 3/2. This confirms that the dynamic exponent z is 3/2. It also confirms the universality with respect 
to the growth rate (0 < s < 1) and the substrate slope (— 1 < m < 1). This is an assuring feature for the KPZ class 
growth models in view of the fact that, in the asymmetric exclusion process with reflecting boundary condition, the 
mass gaps remain finite even in the bulk limit [14]. The discrete Fourier transform of Eq. (83) gives the structure 
factor of finite size system in the scaling form 

S(k, t) = e- iksmT S{s\fl - m 2 |r|fc 3 / 2 , kN) (84) 

for small k as expected from other theoretical grounds [5]. The phase factor e - iksmT in Eq. (84) is due to the imaginary 
term of the mass gaps and is related to the steady state current as discussed in [15] and [5]. Presence of such a phase 
factor is also consistent with the invariance of the KPZ equation under the Galilean transformation which tilts the 
substrate by an infinitesimal m and changes the spatial coordinate by r — ► r — suit [2,25]. Thus, the modulating 
factor in Eq. (83) disappears along the ray r + suit = in space-time and the imaginary part of Eq. (81) offsets this 
shift in space coordinate. 

Henkel and Schiitz [15] showed that the mass gaps of Eq. (1) scale as l/N 2 for q ~ 0(1/ N), i.e. for finite 
number of down steps instead of finite density of down steps. It is consistent with Eq. (81) due to the factor 
\/l — m 2 = y/4q(l — q). Higher order terms in Eq. (81) also contribute to 0(1/N 2 ) when q ~ 0(1/ N). Even if the 
l/N 2 scaling for q ~ 0(1/ N) is the same as that for the s = EW point, their physical origins are different. The 
latter comes from the fact that correlation functions such as Eq. (82) satisfy the discrete diffusion equation exactly 
for all m while the former originates from the non-relativistic free fermion dispersion relation near the band edge 
[3,5,6]. When A < 1, the limits q — ► 0, 1 are where the Pokrovsky-Talapov (PT) type commensurate-incommensurate 
transitions occur in the context of the domain wall picture as discussed in detail in [9] for the case of s = 1. The 
l/N 2 scaling of the mass gap near q = is a signature of the PT transition in this case. 

The second main result of this paper addresses the KPZ to EW crossover scaling function. The mass gaps Eq. (81) 
depend linearly with the growth rate s to the leading orders whereas, when s — ► along the stochastic line, they are 
expected to cross over from the KPZ scaling to the EW scaling. Such non-trivial s-dependence comes from higher 
order terms in the perturbation series. Thus we have pushed the perturbative calculation to 0(1/ N 4 ) for m = (See 
Eqs. (74) and (76)). From the result, the KPZ to EW crossover scaling variable is found to be u = sy/W = sN 2 ~ z 
and the mass gaps are shown to take the scaling form 

E l N = Gi(sN 2 - z )/N 2 (85) 

for m = where the KPZ to EW crossover scaling function Qi(u) behaves as ~ u for u —> oo and ^(0)=constant. 
The perturbative expansion scheme of this work allows one to obtain the expansion of the scaling function Qi(u) in 
powers of 1/u. Our 0(1/ N 4 ) result determines the first three terms in this expansion. The result is shown in Eq. (79) 
for the first excited level. For other levels, only the numerical values of the coefficients in the series change. Derivation 
of the exact form of Qi(u) is an open problem. Also, to what extent the results of this work remain valid in other 
KPZ class models as universal features is an interesting open question. 
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APPENDIX A: RELATION TO ASYMMETRIC SIX VERTEX MODEL 

In Appendix A, we briefly review, for the sake of self- containment, relevant results of the asymmetric six-vertex 
model and derive the asymmetric XXZ hamiltonian as an anisotropic limit of it. Let W{, i=l, . . ., 6, denote the six 
Boltzmann weights of the six- vertex model in the standard labeling [17,26] . We denote by T the row-to-row transfer 
matrix of the general six-vertex model on N x N square lattice under periodic boundary conditions. The number 
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of down-arrows, Q, in each row is conserved and there are four independent parameters parametrizing u> 8 's. The 
eigenvalue A of T in the sector of Q down-arrows is given by 



Q Q 

A N-Q TT/ , W 5 W 6 Zj N _ Q -pr W 5 W 6 . , , 

A = w 1 \\{w 3 -\ — ) + w 4 \\(w 2 ). (Al) 

.7=1 J .1=1 J 



The sets of fugacities {zj} are given by the Bethe ansatz equation Eq. (4) where in the present 



case, 



is the horizontal field and 



H = | \og{wiw 3 /w 2 w A ) (A2) 



2(wiW2W3W4) 1 / 2 



is the standard interaction parameter. It is also convenient to define 



~ _ A _ wiw 2 + W3W4 - w 5 w 6 , A4 . 
cosh(2if) W1W3 + W2W4 



and 



s = tanh(2F) = WlW3 " W2W4 (A5) 
W1KJ3 + W2W4 

for later uses. Note that Zj in Eq. (Al) is the inverse of that defined in [17] and that m = 1 — 2Q / N = 1 — 2q of this 
work corresponds to y in standard six vertex notation. Eigenvectors of T depend only on H and A and hence T form 
a two-parameter family of commuting matrices; i.e. T's with different values of w 8 's commute with each other as long 
as they have the same values of H and A. This is a direct consequence of the fact that the Yang-Baxter equation has a 
solution under the condition of constant A and H [26]. When Wi's assume the limiting values = = = w® = 1 
and u>3 = u>4 = 0, T becomes the shift operator. We call this the anisotropic limit. Suppose Wi are parametrized by 
u and v in addition to H and A in such a way that the limit u = v = corresponds to the anisotropic limit. An 
example of such parametrization for A > 1 is 

wi = exp(i>) sinh(« + v)/ sinh(V ) 
u>2 = exp(— v) sinh(« + v)/ sinh(V ) 
u>3 = exp(2if — v) sinh(w)/ sinh(V) 
W4 = exp(— 2H + v) sinh(«) / sinh(f) 
w 5 = w 6 = 1 

with A = cosh(V) but actual form of parametrization is not important here. 

The quantum chain hamiltonians commuting with T are obtained by taking derivatives of — log T with respect to 
11 or u at u = d = 0. A straightforward calculation yields [26] 

H , = SlogT 



du 

N 



u=v=0 



w\ + m — Wc. — w 



2 5 6 (1 - <7?<t? +1 ) - w' 3 ata- +1 - w' 4 a-at +1 } - w[(N - Q) - w' 2 Q (A6) 



=1 



for the sector Q, where <rf are the Pauli spin operators and w\ stands for the derivative of with respect to u at 
u = v = 0. Corresponding eigenvalue E' of W is, from Eq. (Al), 

, _ d log A 



du 
Q 

J2i( w i + »2"»5- O ~ vo'szj 1 ~ <*j} ~ <(N -Q)- w' 2 Q. (A7) 
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For H! to commute with T, the derivatives should be taken with H and A fixed. This introduces two constraints; 



w[ + w' 2 - w' 5 - w' 6 = A(w' 3 + w' A ) (A8a) 
w' 3 = exp(4H)w' 4 . (A8b) 

Using Eq. (A8) in Eqs. (A6) and neglecting the constant terms, one then easily sees that W becomes the asymmetric 
XXZ hamiltonian Eq. (1) up to a normalization factor (w' 3 + w' A ). Therefore the energy of Eq. (1) is obtained from 
Eqs. (A7) and (A8) up to the same normalization factor and the constants, which is exactly Eq. (3). 

APPENDIX B: PROPERTIES OF Y M 

In Appendix B, we discuss some properties of Y^(y) and Y m (y) defined in Eqs. (40) and (39), respectively. Here, 
Y^ refers to the scaling function associated with the ground state, while Y m refers to that with the first excited level. 
The most useful property is the relation 

YM = ~Y m . 2 (y), (Bl) 

which also holds for other Y m 's associated with other levels including Y^. 

First we consider the case of m odd. Y^ and Y m are real for real y. As y — ► oo, their asymptotic forms are 

m ( ' y U + 2 12t/ 2 + J ' ( ' 



ym ^ ( _ 1) (-)/ V -^(^-^ + ...) ) (B3) 

respectively. To see how Y^(y) behaves as y — ► — oo, we expand Eq. (40) as a series in —l/y paying attention to the 
branch cut of the square root. After a little algebra, we obtain 



&>2,even 

where 



«<»> = £ (-i)./. "*"-')-^-"*') -/.)(-!) 



k-m/2 



(B4) 



To calculate the integral above, we apply the sum formula Eq. (26) to the trivial sum X2/=oO' — l/2) & an d find 
that Jfc = k/(k+ 1) for k even. Thus all coefficients of the series Eq. (B4) are identically zero. This suggests that Y^ 
has an essential singularity at y = — oo and approaches to zero exponentially as y — ► — oo. For Y m (y) we have 

Y m (y) - -2(-,r/ 2 { i _ H^ZII + . . .} (B6) 

as y —> — oo. Numerically, we find that Y° and — Y 3 ° both are positive monotonic increasing function while Yi(y) is a 
monotonic increasing real function passing through at t/o = 1.119066880828 . . .. Due to the relation Eq. (Bl), t/o is 
the position of the maximum of Y 3 (y). 

For m even, one can show that Y£ = for all m using the value of the integral while Y m is a polynomial in y of 
order (m — 2)/2. Specifically, Y- 2 = 2i, Y4 = —Aiy, Ye = i(6y 2 — 2), etc. and Y m ~ im(—y)^ m ~ 2 ^l 2 for y large. 

For other energy levels than the first excited state we are considering here, the last terms of Eq. (39) are different. 
Thus, the numerical value of t/o where Y\ vanishes is different from level to level. For example, Y\ defined in Eq. (80) 
has a zero at t/o = 1.65873919064 . . . However, one can show in general that Y 2 is related to the total momentum P 
of Eq. (14) by 
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FIG. 1. Locus of x = Z^(<f>) for — wq < <f> < wq in the complex x plane for the case of q = 1/2. 



FIG. 2. The scaling function of the first excited energy level. 
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